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PREFACE 


ThW-volUme •  the  first  part  of  a  series  giving  data  for  design  against  acoustic 
fatigue,  has  been  prepared  in  order  to  draw  together  the  results  of  research  in 
acoustic  fatigue  and  to  present  then  in  a  fora  directly  useable  in  aerospaoe  design. 
Future  work  in  this  series  will  deal  with  endurance  of  aluminium  alloy  and  titanium 
alloy  structures  under  simulated  acoustic  loading)  stress  response  of  flat  or  curved 
honeycomb  panels,  near  field  compressor  noise  estimation)  stress  response  of  box  and 
control  surface  structures  and  structural  damping*  ,  j 

The  AOARO  Structures  and  Materials  Panel  has  for  many  years  been  active  in 
encouraging  and  coordinating  the  work  that  has  been  neoessary  to  make  this  oolleotion 
of  design  data  possible  and  after  agreeing  on  procedures  for  the  acquisition)  analysis 
and  interpretation  of  the  requisite  data)  work  on  this  series  of  design  data  sheets 
was  initiated  in  1970. 

The  overall  management  of  the  project  has  been  conducted  by  the  Vorking  Group  >n 
Acoustic  Fatigue  of  the  AGARD  Structures  and  Materials  Panel)  and  the  project  has  been 
financed  through  a  collective  fund  established  by  the  Nations  collaborating  in  the 
project)  namely  Canada)  France,  Oermany)  Italy)  U.X.  and  U.S.  National  Coordinators 
appointed  by  each  oountry  havo  provided  the  basio  data)  liaised  with  the  souroes  of 
the  data,  and  provided  constructive  comment  on  draft  data  sheets.  These  Coordinators 
are  Dr  G.M.  Lindberg  (Canada),  Mr  R.  Loubet  (France),  Mr  0.  Bayerdorfer  (Germany), 
Gen.  A.  Oriselll  (Italy),  Mr  N.A.  Townsend  (U.K.),  Mr  A.V.  Kolb  (U.S.)  and 
Mr  F.F.  Rudder  (U.S.).  Staff  of  the  Engineering  Sciences  Data  Unit  Ltd,  London,  have 
analysed  the  basic  data  and  prepared  and  edited  the  resultant  data  sheets  with  invalu¬ 
able  guidance  and  advice  from  the  National  Coordinators  and  from  tbs  Acoustic  Fatigue 
Panel  of  the  Royal  Aeronautical  Society  which  has  the  following  constitution! 

Professor  B.L.  Clarkson  (Chairman) ,  Mr  D.C.O.  Eaton,  Mr  J.A.  Hay,  Mr  V.T.  Kirkby, 

Mr  M.J.T.  Smith  and  Mr  N.A.  Townsend.  The  members  of  Staff  of  the  Engineering 
Sciences  Data  Unit  concerned  with  the  preparation  of  the  data  sheets  in  this  volume 
arei  Mr  A.G.R.  Thomson  (Executive,  Environmental  Projects),  Dr  G.  Sen  Gupta  and 
Mr  R.F.  Lambert  (Environmental  Projects  Group). 

Data  sheets  based  on  this  AOARDograph  will  subsequently  be  issued  in  the  Fatigue 
Series  of  Engineering  Sciences  Data  issued  by  ESDU  Ltd,  where  additions  and  amendments 
will  be  made  to  maintain  their  current  applicability. 


A.K.  Hall 
Chairman, 

Vorking  Group  on  Acoustic  Fatigue 
Structures  and  Materials  Panel 
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Section  1 

INTRODUCTION  TO  DSSZON  INFORMATION 
ON  ACOUSTIC  FATIGUE 


1*1  Osnsral  Remarks 

Experience  on  many  different  aircraft  haa  shown  that  the  wide  band  high  in tonal ty 
nolee  aaeoclated  with  a  Jet  exhaust  can  cause  structural  fatigue  failure  In  regions 
close  to  the  Jet*  Similar  failures  have  oocurred  In  other  regions  of  pressure  fluctu¬ 
ation  suoh  as  dose  to  propeller  tips  and  in  regions  of  separated  flow*  for  example 
behind  airbrakes*  To  date  there  is  no  known  case  of  catastrophic  failure  of  an 
aircraft  due  to  acoustic  fatigue  but  the  daaage  to  the  structure  can  lead  to  unaccept¬ 
able  expense  for  maintenance*  inspection  and  loes  of  aircraft  use* 

This  Introduction  outlinee  the  important  faotors  in  aoouatio  fatigue  failure  and 
gives  guidance  on  the  use  of  data  in  this  AOARDograph  and  in  relevant  Engineering 
Sciences  Data  Items*  In  the  present  state  of  knowledge  a  complete  analytical  solution 
can  not  be  preeented*  but  the  framework  of  a  design  procedure  applicable  especially 
to  skin  panels  can  be  described*  * 

Terms  used  in  acoustic  fatigue  analysis  are  defined  in  Reference  1*5*2  and  a  short 
bibliography  of  relevant  literature  ia  given  in  Reference  1*5*3* 

1*2  Outline  of  the  Problem 

Table  1*1  summarises  the  important  factors  affecting  the  aoouatio  fatigue  life 
of  a  structure  subjected  to  Jet  noise*  Similar  faotors  apply  for  other  aourcea  of 
pressure  fluctuation*  when  the  relevant  factors  of  importance  for  the  operating 
conditions  and  noise  field  characteristics  need  to  be  considered  in  place  of  1  and  2 
in  Table  1*1* 


Table  l.i 


Principal  Items 

Factors  of  Importance 

• 

a 

o 

■H 

Q 

*0 

fl 

ti 

5 

1*  Operating 
condition 

Flight  plan 

Time  at  full  thrust  and  low  forward  apeed 

Use  of  reverse  thrust 

Use  of  afterburner 

Occurrence  of  shock  cell  noise 

2*  Noise  field 

characteristics 

Engine  characteristics 

Engine  position 

Reflected  noise  from  ground  and  other 
structure 

3*  Structural  response 

Natural  frequencies 

Mode  shapes 

Nodal  damping 

2**  Stress  at  critical 

points 

Predominant  nodes 

Detail  design  (stress  concentration) 

5*  Fatigue  damage 

Long  life  portion  of  S  -N  curve  under 
random  loading 

Crack  propagation  characteristics 

Environment  (oorrosion*  temperature) 

1*3  D»«lwi  froe«dur« 


1.3.1  Skin  panels 

The  basic  type  of  akin  design  for  a  component  aubjaot  to  aoouatlo  fatigue 
loading  la  ohoaan  with  regard  to  the  overall  nolae  level  whioh  may  be  eatimated  aa 
outlined  In  Section  1.4.2  or  from  previous  experience.  Ouidanoe  on  the  oholoe  of  akin 
dealgn  and  on  good  dealgn  practice  is  given  In  Reference  1.5.4.  The  effect  of  variationa 
In  panel  geometry  may  be  lnveatigated  by  means  of  tha  analysis  outlined  In  Seotlons 
1.4*1  to  1.4*4.  Use  of  early  design  charts  auoh  as  those  of  Reference  1.5.1  has  been 
found  to  lead  to  exoesaively  heavy  designs  In  many  oases,  but  the  later  References 
are  more  realistic. 

Toots  of  critical  seotlons  of  strueturs  in  a  simulated  jet  noise  environment  may 
he  neceaaary  after  preliminary  design,  and  after  redesign  a  proof  test  may  be  required 
before  finalisation. 

1.4  Analysis 

The  analytical  estimation  of  the  life  of  a  part  subject  to  acoustic  fatigue 
requires  inf ormc.tion  in  three  main  areas  1 

(a)  the  loading  action,  l.e.  the  nolae  field  cha.  'icteriatics  and  durations 
under  critical  operating  conditions. 

(b)  the  structural  response,  i.e.  the  natural  frequencies  of  the  structure 
and  the  amplitudes  of  the  induced  stresses. 

(c)  the  fatigue  life  of  the  type  of  structure  made  from  the  particular 
material  when  subjected  to  the  resultant  random  loading. 

The  following  notes  amplify  these  information  requirements. 

1.4*1  Operating  conditions 

The  starting  point  for  an  analysis  is  the  estimation  of  the  pressure  fluctuations 
and  their  durations  throughout  the  flight  plan  or  mission  profile,  and  the  estimated 
number  of  flights  during  service  life.  For  example  Important  elements  may  be  the 
total  time  during  which  the  engine  is  at  full  thrust  and  the  forward  speed  is  low.  k 

the  possible  occurrence  of  shock  cell  noise  in  orulse.  and  the  use  of  reveres  thrust 
in  landing.  In  addition  consideration  should  be  given  to  pressure  fluctuations  such 
as  those  generated  by  Jet  exhaust  impingement  and  boundary  layer  turbulenoe. 

Consideration  of  such  elements  leads  to  an  estimate  of  the  acoustic  fatigue  load 
spectrum. 

Reference  1.3.11  contains  a  suggested  procedure  for  dealing  with  cases  where  more 
than  one  operating  condition  is  a  significant  source  of  acoustic  fatigue  damage. 

Reference  1.5.6  deals  with  the  combination  of  sound  levels  from  two  or  more  sources 
acting  simultaneously. 

1.4*2  Noise-field  characteristics 


The  loading  action  on  the  structure  is  computed  from  a  knowledge  of  the  engine 
noise  characteristics  and  tha  vehicle  configuration.  The  most  Important  factors  are 
the  location  of  the  engines  relative  to  the  structure  and  the  proximity  of  sound 
reflecting  surfaces  such  as  wing,  fuselage,  tall  section  or  runwsy.  The  characteristics 
required  to  define  the  loading  action  completely  are 

(a)  the  preaaure  spectra  and  spatial  distribution  over  the  structural 
surface. 

(b)  the  correlation  of  the  pressures  over  the  structural  surface. 

A  method  of  estimating  the  first,  in  relation  to  Jet  noise,  is  giveh  in  Section  2 
and  further  methods  are  listed  in  Reference  1.5.3.  Although  estimates  by  the  method  of 
Section  2  have  compared  well  with  measurements  on  several  different  full-scale  and 
model  aircraft,  it  usually  remaino  necessary  in  the  design  development  stage  to 
supplement  them  by  measurements  of  noise  levels  on  prototype  aircraft. 

The  spatial  correlation  of  the  pressures  is  not  usually  required  in  initial 
design  because  of  the  simplifying  assumptions  which  must  be  made. 

1.40  Structural  response  and  stress  at  critical  points 

The  complete  expression  for  the  response  of  a  structure  to  acoustic  excitation 
ruquires  n  knowledge  of  all  the  modes  of  vibration  of  the  struoture.  their  natural 
frequencies  and  damping,  in  addition  to  the  noise  field  characteristics.  In  praotioe. 
the  response  of  the  structure  can  be  estimated  only  from  simplified  representations  of 
the  struoture  and  of  the  pressure  field. 
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The  following  table  ahowo  sources  of  data  on  natural  frequencies  of  various  typos 
of  structure  in  their  lowor  modes • 


Type  of  Structure 

Edge  Conditions 

Data  Source 

Plot  piote 

All  edges  fixed 

All  edges  simply-supported 

Ref.  1.5.8 

O 

Flot  honeycoinb  pone) 

All  edges  simply-supported 

Ref  1  5.9 

(to  be  revised  in  1972) 

Stringer-stiffened  ponri 

Torc-ionol  restraint  ot  inflexible  stringers, 

long  edge  of  skin  simply-supported, 

short  edge  restrained  in  rototion 

Section  3. 

Curved  plate 

All  edges  fixed 

All  edges  simply -supported 

(0  >  b) 

Section  4 

Curved  honeycomb  ponet 

All  edges  fixed 

All  edges  simply -supported 

To  be  issued  in  Port  2, 1972 

Zn  calculating  natural  frequencies  and  stress  response,  experience  suggests  that 
individual  plate  widths  should  be  taken  as  tho  distsnce  between  adjacent  skin-stiffener 
rivet  lines.  For  bonded  structures,  the  distance  between  centrelines  of  adjacent 
bonded  flanges  is  genarally  taken. 

Knowing  the  natural  frequencies,  tho  Bound  pressure  levels  and  the  pressure 
spectra,  and  making  assumptions  on  panel  damping.  Section  5  gives  an  approximate 
mothod  of  estimating  the  stress  response,  in  terms  of  r.m,s,  stress,  of  skin-stringer 
panels  having  stiffeners  of  relatively  high  flexural  stiffness;  Reference  1.12 
outlines  the  simplified  theory  used  and  shows  a  comparison  with  experimental  results. 

In  using  this  mothod  it  is  required  to  know  the  spectrum  level  of  acoustiu  pressure  at 
a  natural  frequency  of  the  panel.  If  band  levels  only  are  known,  they  may  be  converted 
to  spectrum  level  using  Reference  1.5«5--  Tho  effect  on  the  stresses  of  reinforcement  by 
increased  thickness  along  the  plate  edges  is  considered  m  Reference  1,5.13. 

1.4»4  Fatigue  strength  and  endurance 

The  estimated  r.m.s.  stress  in  the  panel  is  compared  with  appropriate  latigue 
strength  data  in  order  to  estimate  the  life  of  the  panel.  In  making  this  comparison 
it  is  important  to  ensure  that  the  fatigue  strength  data  used  correspond  to  the 
appropriate  conditions  because  the  highest  stresses  usually  occur  in  tne  region  of 
stiffener  or  frame  rivet  lines  where  effects  of  stress  concentration  and  fretting  are 
important.  Reference  1.10  gives  data  on  the  endurance  of  var >oup  types  of  Joint  in 
aluminium  alloys  subjected  to  simulated  acoustic  loading. 

Where  a  panel  is  subjected  to  acoustic  loadings  of  different  intensities  at 

different  times  it  is  suggested  tentatively  that  the  effects  should  be  combined  by 

applying  the  simple  cumulative  damage  hypothesis  using  S  -N  data.  This  approach 

roB  r 

is  discussed  in  Reference  1.5,11. 

1.5  References 
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Engineering  Sciences  Data  Item  Ho.  66013,  1966. 
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1.5.11  Kirkby,  V.T. 

Edwards.  P.R. 


1.5.12 
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1.5.14 
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Engineering  Sciences  Data  Item  No.  66014,  1966. 

General  principles  for  the  doaign  of  acoustically 
or  similarly  excited  structure. 

Engineering  Sciences  Data  Item  No.  66015>  1966. 

Bandwidth  correction. 

Engineering  Sciences  Data  Item  No.  66016,  1 966. 
Combination  of  levels  in  dB. 

Engineering  Sciences  Data  Item  No.  660!7»  1966. 

The  relation  between  sound  pressure  level  and 
r.m.s.  fluctuating  pressure. 

Engineering  Sciences  Data  Item  No.  66018.  1966. 

Natural  frequencies  of  uniform  flat  plates. 

Engineering  Sciences  Data  Item  No.  66019.  1966. 

Natural  frequencies  of  flat  sandwich  panels  with 
cores  of  zero  flexural  stiffness  and  simply- 
supported  edges. 

Engineering  Sciences  Data  Item  No.  66020.  1966. 

Endurance  of  riveted  skin-rib  flange  connections 
(tentative).  (Aluminium  alloy  material  -  In  bonding.) 
Engineering  Sciences  Data  Item  No.  66022,  1966. 

(To  be  superseded  in  Part  2  of  this  AGARDograph.) 

A  method  of  fatigue  life  prediction  using  data 
obtained  under  random  loading  conditio  s. 

RAE  tech.  Rep.  66023.  1966. 

Estimation  of  the  r.m.s.  stress  in  skin  panels 
subjected  to  random  acoustic  loading. 

Engineering  Sciences  Data  Itom  No.  67028,  1967* 

The  effect  of  edge  reinforcomont  on  the  strecseo 
in  skin  panels  under  uniform  pressure. 

Engineering  Scionces  Data  Itom  No.  67029,  1967. 

Hear-field  noise  analysis  of  aircraft  propulsion 
systems  with  emphasis  on  prediction  techniques 
for  Jets. 

AFFDL-TR-67-43*  1967. 

Refinement  of  sonic  fatigue  structural  design  criteria. 
AFFDL-TR-67-156,  1968. 


5 


Section  2 

THE  ESTIMATION  OP  NEAR  FIELD  SOUND  PRESSURE 
LEVELS  DUE  TO  JET  NOISE 


2.1  Notation 


°p 

specific  heat  at  constant  pressure  of 
fully  expanded  jet  gases 

j/kg  K 

ft  lbf/alug  K 

D 

Jet  nozzle  diameter 

a 

ft 

f 

frequency 

Hz 

c/e 

L 

overall  sound  pressuro  level 

dB 

dB 

Lo 

datum  overall  sound  pressure  level 

dB 

dB 

Lv 

correction  to  L  due  to  change  m  jet 
velocity 

dB 

dB 

Lp 

correction  to  L  due  to  change  in  jet 
density 

dB 

dB 

n 

velocity  index 

P 

ambient  pressure 

N/m2 

lbf/ft2 

R 

gas  constant 

(287  J/kg  K)  (3090  ft  lbf/alug  K) 

J/kg  K 

ft  lbf/slug  K 

TJ 

Jet  pipe  static  temperature 

K 

K 

thj 

Jet  pipe  total  temperature 

K 

K 

V 

mean  fully  expanded  Jet  velocity 

m/o 

ft/s 

Vo 

datum  fully  expanded  Jet  velocity 
(610  m/e)  (2000  ft/s) 

ra/  s 

ft/s 

X 

axial  distance  from  jet  pipe  nozzle 
plane,  measured  in  jot  direction 

m 

ft 

y 

radial  diatance  from  Jet  axis 

D 

ft 

pJ 

density  of  Jet  gases 

kg/m3 

9 lug/ft3 

Po 

datum  density 

(0.A9  kg/a?)  (0.00C  35  slug/ft3) 

kg/m3 

slug/ft3 

Both 

SI  and  British  units  are  quoted  but  any  coherent 

system  of 

units  may  bo  used. 

2.2  Introduction 

This  Section  gives  a  method  of  estimating  the  near  field  sound  pressure  le\els 
due  to  high  »elfrcity  Jet  noise  from  a  single  stationary  conical  nozzle. 

Detailed  noise  data  from  the  engine  manufacturer  should  be  used  for  preference 
but,  when  this  is  not  possible,  a  number  of  methods  of  estimation  is  available.  These 
methods  giv*  widely  differing  results  and  after  comparing  them  the  following  method 
has  been  found  Co  give  the  closest  agreement  with  experimental  results.  This  method 
is  intended  to  give  a  value  of  the  overall  sound  pressure  level  that  is  somewhere  near 
the  mean.  Tht  accuracy  of  tue  estimated  V-vels  is  expected  to  be  within  -5  dB,  ) 
betvoen  approximately  60  per  sent  and  17’-  per  cent  of  the  actual  r.m.s.  preasure. 

Under  reheat  condition'*  or  at  y  i  ^elocitd  *t>  woll  above  the  datum  velocity,  however, 
the  discrepancy  emit'  he  .Luvri 

Tho  noiee  lavel*  prsdic* .d  by  the  method  described  here  &r  free  field  values  but 
do  not  apply  to  the  field  within  a  conical  ourface  emanating  fr«,  the  jet  nozzle  peri¬ 
meter  and  expanding  at  a  cemi-angJe  of  approximately  15  degrees  i  lative  to  the  jot 
axis,  since  this  region  contains  the  turbulent  mixing  region  of  tue  jet. 

This  Section  is  also  inapplicable  to  the  case  of  "shock  col}  no  is*. " ,  u 
phenomenon  ufaich  oan  occur  when  a  jet  nozzle  is  operating  in  an  over -choked  condition. 
A  relatively  discrete  note  is  then  produced,  giving  in  a  nurrow  band  an  intensity  wuct 
greater  than  that  of  the  basic  jet  noise. 
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2.3  Calculation  Procedure 

The  following  information  is  required  a <t  a  basis  for  the  computation  of  the  sound 
pressure  level  at  seme  point  in  the  field! 

jet  nozzle  diameter 
co-ordinates  of  the  point 
mean  fully  expandod  Jet  velocity 
jet  pipe  total  temperature 

(i)  Evaluate  x/D  and  y/D  • 

(ii)  From  Figure  1,  which  shows  the  free  field  datum  noise  level  contours , 
read  off  the  datum  overall  sound  pressure  level  L  at  the  position 
(x/D,  y/D). 

(iii)  From  Figure  2,  which  shows  volocity  index  contours,  read  off  the  velocity 
index  n  at  the  position  (x/D,  y/D). 

(iv)  Calculate  the  first  part  of  the  velocity  correction  term  ALyl  ,  from 
ALvl  -  lOn  log1Q  . 

(v)  Obtain  tho  second  part  of  the  velocity  correction  term  AL  2  from 

Figure  3,  which  shows  AL  2  plotted  against  log1Q(V/v7)  V  for  different 
values  of  n  •  J 

(vi)  Calculate  the  total  jet  velocity  correction  Ly  from 

Lv  -  Alv1  ♦  Alv2  . 

(vii)  Calculato  the  density  of  the  jet  gases  Pj  , 
whoro  Pj  .  -E_ 

and  T j  «>  THj  -  for  all  jet  velocities. 

(viii)  Calculate  donsit/  correction  Lp  , 

where  Lp  -  801o*xo(^)* 

(ix)  Evalunto  tho  overall  sound  pressure  level  (SPL)  at  »  required  point 
(x,  y)  as  * 

L  w  1  ♦  L  +  L  . 

o  v  p 

This  level  can  be  converted  to  a  pressure  loading  by  the  use  of  the 
relation 

r.m.s.  fluctuating  pressure  in  N/m2  »  j0(0.05L-4.699) 
or  using  Reference  2.3.1. 

(x)  For  a  given  frequency  the  spectrum  level  in  decibels  relative  to  an 

arbitrary  overall  SPL  may  be  read  from  Figure  It,  after  evaluating  fD/V 
and  x/D  .  Alternatively,  if  the  spectrum  level  is  required  in  terms  of 
pressure  loading,  this  may  be  found  using  Figure  5*  Taking  values  of 

fD/V  and  x/D  the  spectrum  level  in  (N/m2)/Hz  relative  to  1  N/m2 
overall  SPL  is  read  froia  the  ordinate  of  Figure  5  and  is  then  multiplied 
by  the  value  of  L  obtained  as  a  pressure  loading  in  (ix).  This  gives 

the  pressure  spectrum  level  for  a  specified  frequency  in  (N/m  )/Hz. 


2,  It  Use  of  Calculated  Sound  Pressure  Levels 


After  the  free  field  sound  pressure  levels  have  been  'etimated,  corrections  must 
be  made  for  local  effects,  for  example  reflection.  When  tho  wavo  fionts  of  this  noise 
field  strike  a  structure  they  are  partially  reflected,  and  the  reflection  process 
locally  increases  the  pressure  loading  on  the  surface. 

When  the  wave  front  strikes  the  surface  at  right  angles  (normal  incidence)  vhe 
pressure  loading  is  doubled  (a  6  dB  Increase  on  the  cumulated  values).  Zf  the  wave 
front  moves  parallel  to  the  surface  (gracing  incidence)  there  is  no  increase  in  loading. 


To  allow  for  this  reflection  process ,  it  is  fairly  general  engineering  practice 
to  add  a  mean  correction  of  3  dB  to  the  calculated  free  air  levels* 

The  effects  of  different  nozzle  arrangements,  aircraft  configurations,  etc., 
should  be  taken  into  consideration. 


Further  corrections  must  be  made  for  the  case  of  an  eircraft  in  motion  to  allow 
for  the  effects  of  the  aircraft  velocity  and  ambient  speed  of  sound. 

An  alternative  computer-based  method  that  plots  overall  sound  pressure  levels  for 
various  engine  operating  conditions  can  be  seen  in  Derivation  2.5*3. 

2.5  Derivation  and  Reference 

Derivation 


2.5*1  Franken,  P.A. 
et  al. 


Methods  of  flight  vehicle  noise  prediction. 
VADC  TR  58-343,  1958. 


2.5.2 


Unpublished  work  by  Rolls-Royce  Ltd  and  Bristol  Siddeley 
Engines  Ltd. 


2.5.3  Plumblee,  H.E. 
et  al. 


Near  field  noise  analyses  of  aircraft  propulsion  systems 
with  emphasis  on  prediction  techniques  for  jots. 

AFFDL-TR-67-43,  1967. 


2.5.4  Riley,  M.P.  Near  field  Jet  noise  prediction  techniques. 

British  Aircraft  Corporation  Ltd,  Acoustics  Laboratory 
Report  A.R.324,  1971. 


Roference 

2.5.5  The  relation  between  sound  pressure  level  and  r.m.s. 

fluctuating  pressure. 

Engineering  Sciences  Data  Item  No.  66018,  1966. 


2.6  Example 


It  is  required  to  estimate  the  r.m.s.  sound  pressure  level  and  the  pressure 
spectrum  level  at  300  Hz  at  a  specific  point  given  the  following  conditions: 


D  ■  0.6l  m, 

V  a  670  m/n, 

V  =  610  m/s, 
o 

R  -  28?  J/kg  K, 
p  n  101  x  103  N/m2. 


x  «  5*5  m» 


y  b  6.7  m, 


Thj  -  900  K, 


p  a  O.49  kg/m3, 


cp  -  1160  J/kg  K, 


Firstly 


x 

D 


0.61 


9.02  , 


and 

Hence,  from  Figure  1, 


Z 

D 


0761 


140  dB 


11.0  . 


and,  from  Figure  2,  by  interpolation 

n  b  6.2  . 


Now 


ALyl  »  10  x  6.2  log10. 


2.5  dB 


and,  from  Figure  3,  for  ^06^q(^/Vo)  «  0.0407  and  interpolating  for  n  *  6.2  , 


so  that 


AL 


'v2  “ 

-  0.9 

dB 

Lv  0 

2,5  - 

0.9  a  1.6 

dB  . 

900  - 

6702 

707  K 

TJ  “ 

2  x  1160  B 

As 


8 


th*"1 

PJ 

* 

HH?§7  *  °*W8  k«/°3  * 

Now 

S 

m 

!01°«10  (^)-°-»  dB 

and  so.  finally, 

L 

» 

140  4  1.6  4  0.14  »  142  dB  . 

To  determine 
Figure  5  ia  uaed. 

the  corresponding  pressure  spectrum  level?  for  a  frequency  of  300  Hz, 

fD 

V 

m 

^70-'-  "  0.273 

and 

X 

D 

- 

9.02  . 

By  into'  olation,  the  spectrum  level  relative  to  an  overall  sound  pressure  level  of 
1  U/n  is  found  to  be  0.02 6  (N/m2)/Hz  • 

2 

In  this  Example  L  *  142  dB  ,  which  from  Paragraph  3( ix)  is  equivalent  to  252  N/m  • 
So  the  pressure  spectrum  level  »  0.026  x  252  «  6.55  (N/m  )/Hz  • 
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Section  3 

NATURAL  FREQUENCIES  OF  BUILT-UP,  FLAT, 
PERIODIC  SKIN  STRINGER  STRUCTURES 

(Part  It  Stringers  rigid  in  bending) 

3*1  Notation 


& 

width  of  akin  plate  or  frame  pitch 

a 

in 

b 

stringer  pitoh 

m 

in 

D 

flexural  rigidity  per  unit  width  of  the  skin 

plate,  Et3/[l2(l-ff2)] 

N  a 

lbf  in 

E 

Young's  modulus  of  skin  material 

N/m2 

lbf/ in' 

Es 

Young's  modulus  of  stringer  material 

N/m2 

lbf/ in‘ 

f 

frequency 

Hz 

c/s 

°s 

shear  modulus  of  stringer  material 

N/m2 

lbf/  in‘ 

polar  moment  of  inertia  of  stringer  cross 
section  about  point  on  skin  directly 
beneath  shear  centre  of  stringer 

m4 

in4 

St  Venant  constant  of  uniform  torsion  for 
stringer  cross  section 

m4 

in4 

K 

frequency  parameter 

m/s 

in/s 

M1 

mode  parameter  for  first  group  of  natural 
frequencies 

m2 

mode  parameter  for  second  group  of  natural 
frequencies 

m 

number  of  half  waves  across  frame  pitch 

N 

nunber  of  spans 

r 

mode  number 

t 

skin  thickness 

m 

in 

V 

velocity  parameter  for  skin  material* 

r. 

warping  constant  of  stringer  cross  section 
with  respect  to  the  point  on  skin  directly 
beneath  shear  centre  of  stringer 

»6 

in6 

*R 

non-dimensional  torsional  stiffness  of 
stringer  »  actual  rotational  stiffness 
per  unit  length  of  the  stringer  x  b/D 

P 

density  of  skin  material 

kg/m3 

K 

P. 

density  of  stringer  material 

kg/m3 

II 

9 

Poisson's  ratio  of  the  skin  material 

Both 

1 

SI  and  British  units  are  quoted  but  any  coherent 

system  of  units 

may  be 

*  The  velocity  parameter  is  defined  in  SI  units  ae  V  ■  (E/p) */2/5080  and  in  British 

unite  as  V  «  (E/p)l^?/200  000,  V  is  approximately  unity  for  all  common  structural 
metallic  materials, 

*  A  density  value  expressed  in  British  units  as  pounds  per  cubic  inch  has  to  be 
divided  by  386, 4  before  it  can  be  used  in  tbs  formula  given  here,  (A  force  of  1  lbf 

2 

acting  on  a  mass  of  1  lb  produces  an  acceleration  of  386,4  in/s  ,) 
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3.2  Introduction 

This  Section  givos  the  natural  frequencies  of  vibration  of  built-up,  flat  skin- 
stringer  structures  with  equispaced  and  identical  stringers  as  shown  in  Figures  3»1 
and  3*2.  The  natural  frequencies  of  such  multi-span  structures  occur  in  groups,  the 
number  of  frequencies  in  each  group  being  equal  to  the  number  of  spans.  For  a  panel 
having  a  large  number  of  spans  the  lowest  natural  frequency  in  any  group  is  generally 
associated  with  a  mode  in  which  the  stringers  twist  and  the  highest  nttural  frequency 
in  that  group  is  associated  with  a  inode  in  which  the  stringers  tend  to  bend  without 
twisting.  For  a  panel  having  a  small  number  of  spans,  typical  mode  shapes,  as  shown 
in  Figure  3*3,  ore  affected  by  the  edge  conditions.  The  mode  shapes  alternate  between 
symmetric  and  anti-symmotric  as  frequency  increases. 

Xn  the  Derivation  the  frame  torsional  stiffness  has  bean  neglected  and  the  frame 
bending  stiffness  has  been  assumed  to  be  infinite.  The  skin  edges  at  the  frames  are 
thus  assumed  to  be  simply-supported.  The  natural  frequencies  are  only  slightly 
influenced  by  the  stringer  bending  stiffness  and  therefore  this  stiffness  bus  been 
assumed  to  be  infinity. 

This  Section  can  be  used  to  predict  the  first  two  groups  of  natural  frequencies 
of  skin-stringer  structures  with  four  different  end  conditions.  The  results  have  been 
presented  for  varying  degrees  of  stringer  torsional  stiffness,  for  three  different 
aspect  ratios  (a/b)  and  for  two  different  values  of  a  ,  the  number  of  half-waves 
across  the  frame  pitch. 

3»3  Calculation  Procedure 

Calculate  f  from 

f  ■  VK  -4  . 

b* 

To  determine  K 

(i)  Calculate  *R  from  the  equation 

"t  ■  {  E.r.  (t?)  4  •  V.(¥)S  -  >.«.<•«*>*} 

Since  f  is  unknown  its  value  has  to  be  estimated  initially.  Values  of  f  *  100  Hz 
for  tho  first  group  of  frequencies  and  f  «  200  Hz  for  the  second  group  should  give 
results  of  sufficient  accuracy  for  commonly  used  skin-stringer  structures.  Alternatively, 
for  the  first  group  of  natural  frequencies,  f  may  be  taken  as  the  fundamental 
natural  frequoncy  of  any  individual  panel  with  fully-fixed  edges.  For  the  second 
group,  f  may  be  taken  as  the  natural  frequency  of  the  same  panel  with  fully-fixed 
edges  in  its  second  mode  (i.e.  with  approximately  two  half  waves  along  the  length  b) . 
Both  these  frequencies  can  be  calculated  using  Reference  3«5<3« 

(ii)  Calculate  Mj  or  Mg  from  the  following  tablet 


(iii)  Read  K  from  Figures  4-9  for  appropriate  values  of  m  ,  aspect  ratio 
(a/b),  kr  and  Mj  or  Hg  . 

Note t  For  structures  with  end  conditions  as  in  Figure  2d,  K  is  approximately 
equal  to  the  arithmetic  average  of  the  values  of  K  obtained  from  the  first  two  rows 
in  the  above  table. 
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3*4  Note* 


The  effect  of  including  the  frame  torsional  stiffnese  would  be  to  Increase  the 
natural  frequencies  slightly  and  this  would  be  more  pronounced  for  structures  with 
smaller  aspect  ratios.  The  effect  of  including  the  framo  bending  stiffness  would  be 
to  decrease  the  natural  frequencies  slightly. 

For  structures  with  a  low  stringer  torsional  etiffnessi  the  bounding  froquencier- 
of  the  first  group  (with  m  *  1  or  2)  can  also  be  found  from  Data  Item  No*  6601?* 

Derivations  3*5*2,  3*5*4  and  3*5*5  give  the  transfer-matrix  method  of  analysis  of 
skin-stringer  structures  which  could  be  used  to  calculate  the  natural  frequencies  of 
non-periodic  structures,  i.e>  structures  with  non-uniform  skin/stringer  characteristics* 

3*5  Derivation 


3*5*1 

Lin,  Y.K. 
Brown,  I.D* 
Deutschle,  P.C 

Free  vibration  of  a  finite  row  of  continuous  skin-stringer 
panels.  J.  Sound  and  Vibration,  Vol.l,  No.l,  January  1964* 

• 

3*5*2 

Lin,  Y.K. 
et  al. 

Free  vibration  of  continuous  skin-stringer  panels  with 
non-uniform  stringer  spacing  and  panel  thickness. 
AFML-TR-64-347, . February  1965* 

3*5*3 

Natural  frequencies  of  uniform  flat  plates. 

Engineering  Sciences  Data  Item  No.  66019,  1966. 

3*5*4 

Mercer,  C.A. 
Seavey,  C. 

Prediction  of  natural  frequencies  and  normal  modes  of 
skin-etringor  panel  rows. 

J.  Sound  and  Vibration,  Vol.6,  No.l,  January  1967. 

3*5-5 

Mercer,  C.A* 
Seavey,  C. 

Program  for  the  calculation  of  natural  frequencies  and 
normal  modes  of  skin-stringer  panel  arrays. 

Institute  of  Sound  and  Vibration  Research  Tech.  Report 

No. 6,  July  1968. 

3*5*6 

Sen  Gupta,  0. 

Natural  flexural  waves  and  the  normal  modes  of 
periodically  supported  beams  and  plates. 

J.  Sound  and  Vibration,  Vol.13,  No.l,  September  1970* 

3*5*7 

Sen  Gupta*  G. 

Natural  frequencies  of  periodic  skin  stringer  structures 

us..ug  a  wave  approach* 

J.  Sound  and  Vibration,  Vol.l6,  No. 4,  June  1971* 


3*6  Example 

It  is  required  to  determine  the  first  group  of  natural  frequencies  of  a  six-span 
akin-etringer  structure  having  the  following  dimensions  >^tc.  with  end  conditions  he 
in  Figure  2d  and  with  a  single  half  wave  across  the  skin  widtn* 

r  ■  3*1  x  10"12  m6 

O 

I8  -  10*5  x  10“8  a* 

E  -  Ef  «  72  300  MN/m2 
p  u  p8  »  2800  kg/nf* 

0.  ■  27  800  MN/m2. 

0 

N  m,  V  »  -  1.00 

5080 

and  for  the  first  group  of  natural  frequencies  with  a  single  half  wave  across  the  skin 
width  (m  m  1  and  assuming  f  a  100  Hz) , 

*r  -  {v»  (!) 4  ♦  Vs  (D  2  -  p.1.  '*  *  100>2 }  * ! 

a  (328  ♦  100  -  116)  x  0.0384 

a  12.0, 

For  this  group  of  frequencies,  with  m  ■  1.  a/b  a  2.0,  ■  12.0  and  using 

Figure  5a,  the  results  can  be  calculated  and  presented  in  the  form  of  the  following 


a  a  508  mm 

b  a  254  tsm 

t  a  1  mm 
tr  a  0.3 

J,  a  94.4  x  10"12  m4 

m 


From  this  D  ■ 


£t- 


2\ 


6.62 


Mode 

nuaber 

End  conditions 
as  in  Figure  2a 

End  conditions 
as  in  Figure  2b 

End  conditions 
as  in  Figure  2d 

f  <*  VK  -4 

b* 

T 

M1 

K 

(from  Figure  5a) 
(m/s) 

M1 

K 

(from  Figure  5*) 
(m/a) 

K 

(Average  of  the 
results  in 
columns  III  and  V) 
(a/s) 

(using  K  from 
column  VI) 

(Hz) 

(1) 

(II) 

(HI) 

(IV) 

(V) 

(VI) 

(VII) 

1 

0 

4270 

1/6 

4360 

4315 

66.9 

2 

1/6 

4360 

1/3 

4610 

4485 

69.5 

3 

1/3 

'*610 

1/2 

4970 

4790 

74.2 

4 

1/2 

1*970 

2/3 

5370 

3170 

60.1 

5 

2/3 

5370 

3/6 

3690 

5530 

85.7 

6 

5/6 

5690 

1 

5330 

5760 

89.3 

Similar  calculations  can  ba  dono  for  a  structure  with  end  conditions  as  in 
Figure  2c.  It  nay  bo  noted  that  the  frequencies  obtained  in  column  VII  are  less  than 
100  Ha  and  therefore  the  value  of  has  been  slightly  under-estimated.  This  means 

the  actual  frequencies  are  slightly  higher  in  this  particular  problem.  Further 
calculation!!  show  that  the  calculated  frequencies  are  within  about  3.5  per  cent  of  the 
frequencies  obtained  by  taking  a  further  iteration. 
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FIGURE  31  A  PERIODIC  SKIN  STRINGER  STRUCTURE 


Stringer 

tortionol 

stiffness 


b. 


Stringer 

torsionol 

stiffness 


e. 


Stringer 
torsiona 
stiffness 

*»  %  **  Kn  *»  «■ 

d. 


mum  ev.ujsTTjajrsrRsreiiCni’cS  WITH  DIFFERENT  BOUNDARY  CONDITIONS 
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FIGURE  3.6.  f-*2-0,  m=l,  GROUP  I 


FIGURE  3.7.  |  =  20,  m*2,GR0UPI 


FIGURE  3.10.  f  r|.5,  m=|,  GROUP  2 


M, 

FIGURE  3.11.  §  =1*5,  m  =  2,  GROUP  2 
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Section  4 

NATURAL  FREQUENCIES  OF  RECTANGULAR 
SINGLY-CURVED  PLATES 


4.1  Notation 


a 

length  of  longer  plate  side 

» 

in 

b 

arc  length  of  shorter  plats  side 

Q 

in 

E 

Young's  modulus  of  plate  material 

N/m2 

lbf/ in' 

f 

natural  frequency  of  vibration  of  plate 

II* 

c/s 

K 

natural  frequency  parameter  for  fixed  edge  plates 

ra/a 

in/e 

Km,n 

natural  frequency  parameter  of  simply-supported 
plate  in  (m,n)th  mode 

m/a 

in/s 

m 

number  of  half  waves  in  plate  parallel  to  longer 
side 

n 

number  of  half  waves  in  plate  parallel  to 
shorter  (curved)  side 

R 

radius  of  curvature  of  plate 

m 

in 

t 

thickness  of  plate 

m 

in 

V 

volocity  parameter  for  plate  material* 

a 

Poisson's  ratio  of  plate  material 

P 

density  of  plate  muterial 

kg/«a3 

N 

Both 

SI  and  British  units  are  quoted  but  any  coherent  system 

of  units 

may  be 

uaed. 


4*2  Notoa 


4.2.1  Platen  with  simply-supported  edges 


The  natural  frequency  is  given  by 

f  -  V  K. 


n2t 
m,n  ^2 


In  Figure  4*1  values  of  K  are  plotted  against  an/mb  for  a  range  of  valuea 

10)  n 

of  b^n^t  for  initially  unstressed,  cylindrically  curved,  rectangular  platea  of 
uniform  thickneaa  having  all  edgea  aimply  supported.  For  plates  having  all  edges 
simply  supported  the  modes  are  exactly  defined  by  the  number  of  half  waves  across 
the  length  and  breadth  of  the  panel.  These  natural  frequencies  can  be  presented  on 
one  figure  since  the  condition  at  a  nodal  line  is  the  samo  as  that  at  the  boundaries 
parallel  to  the  nodal  line. 

4.2.2  Plates  with  fixed  edges 

The  natural  frequency  is  given  by 

f  -  V  K  ~  • 

b* 


*  The  velocity  parameter  is  defined  in  SI  units  by  V  «  (E/p)fy5080  0r  in  British 

units  by  V  »  (E/p)fy200  000.  V  is  approximately  unity  for  all  common  structural 
metallic  materials. 

*  A  density  value  expressed  in  British  units  at  pounds  per  cubic  inch  has  to  be 
divided  by  386.4  before  it  can  be  used  in  the  formula  for  V  given  here.  (A  force 

of  1  lbf  acting  on  a  mass  of  1  lb  produces  an  acceleration  of  386.4  in/s  . ) 
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In  Figures  4«2  to  4*8  values  of  K  are  plotted  against  a/b  for  a  range  of 

o 

values  of  b  /Rt  for  a  number  of  natural  modes  for  initially  unstressed,  cyxlndrioally 
curved,  rectangular  plates  of  uniform  thickness  having  all  edges  fully  fixed.  These 
modes  of  vibration  are  not  necessarily  uniquely  identified  by  the  number  of  half  waves 
across  the  length  and  breadth  of  the  plate.  In  this  Section  the  fixed-edge  plate 
modes  are  defined  by  the  symmetry  or  antisymmetry  in  two  directions!  stating  first 
the  symmetry  condition  parallel  to  the  plate  longer  side  followed  by  the  symmetry 
condition  parallel  to  the  plate  shorter  curved  side.  The  order  of  modes  in  any  mode 
type  is  given  in  ascending  order  of  frequency. 

The  fixed  edge  panel  modes  condidered  in  this  Section  are  listed  in  Table  I. 


Table  I 


Mode 

Figure  for 

frequency 

parameter 

Figure  for 
mode  shape 

First  (S,S) 

4.2 

4.9 

Second  (S,S) 

4.3 

4.10 

First  (S,A) 

4*4 

4.11 

Second  (S,A) 

4.5 

4.12 

First  (A,S) 

4*6 

4.13 

Second  (A,S) 

4.7 

4.14 

First  (A, A) 

4.8 

4.15 

S  n  symmetric  A  ■  antisymmetric 


The  expression  used  to  calculate  the  fixed-edge  plate  frequencies  contains 
parameters  a/b,  ba/Rt  and  b/R  .  The  curves  in  this  section  are  for  a  value  of 
b/R  of  0.1,  but  tho  curves  may  be  used  for  shallow  plates  generally.  At  a  value  of 
b/R  of  1.0  the  frequencies  are  2-3  per  oent  less  than  obtained  from  this  Section. 

4*2.3  General  notes 


In  deriving  the  curves  the  value  of  O’  was  assumed  to  be  0«3»  This  value  is 
sufficiently  accurate  for  all  common  structural  metallic  materials. 

In  the  theory  on  which  this  Section  is  based  it  is  assumed  that  panels  are  thin 
and  of  shallow  curvature.  These  conditions  are  satisfied  when  R/t>20  and  b/R  <1.5* 

Shear  deflection  and  rotary  inertia  have  been  neglected  when  calculating  the 
curves  for  both  simply-supported  and  fixed-edge  plates.  This  is  a  generally  accepted 
assumption  if  the  ratio  of  half  wave  lengths  to  plate  thickness  is  greater  than  10. 

For  flat  plates  with  simply-supported  edges  the  lowest  natural  frequency  of 
vibration  occurs  in  the  mode  having  a  single  half  wave  across  the  length  and  breadth 
of  the  plate,  i.e.  the  (l,l)  mode.  For  flat  plates  with  fixed  edges  the  lowest 
natural  frequency  occurs  in  the  first  symmetric  -  symmetric  mode  where  the  mode  shape 
is  similar  to  the  simply-supported  edge  (l,l)  mode.  As  the  curvature  of  the  plate 
increases  (radius  reduces)  the  lowest  frequency,  for  both  simply-supported  and  fixed- 
edge  plates,  tends  to  occur  in  a  mode  where  the  number  of  half  waves  on  the  curved 
side  is  greater  than  one. 

4 .3  Application  to  Aircraft  Structural  Panels 

In  typical  aircraft  construction,  stiffening  members  such  as  ribs,  frames  and 
stringers  divide  the  panel  into  plates.  Owing  to  the  mechanical  coupling  between 
plates,  normal  modes  result  in  which  arrays  of  plates  vibrate  together.  Theoretical 
analysis  Indicates  that  an  array  of  plates  has  groups  of  resonant  frequencies.  Each 
group  is  bounded  by  a  lower  frequency  mode  where  stringers  twist  and  a  higher 
frequency  mode  where  stringers  bend.  These  lower  and  higher  frequency  modes  correspond 
to  the  single  plate  modes  with  simply-supported  and  fixed  edges  respectively.  In 
between  these  two  bounding  frequencies  there  is  a  number  of  intermediate  modes  in 
which  the  stringers  both  twist  and  bend.  The  second  and  subsequent  frequenoy  groups 
correspond  to  overtone  modes  of  the  fundamental  group. 

The  number  of  plates  showing  predominant  vibration  at  any  one  natural  frequenoy 
of  the  group  depends  on  the  relative  dimensi.^.  -f  the  plates  and  the  phase  relation¬ 
ships  of  the  exciting  force  over  the  surface  of  the  panel. 


i 


„l 
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The  lowest  natural  frequency  le  generally  associated  with  a  mode  of  vibration 
in  which  the  response  of  plates  is  such  that  half  waves  on  the  two  sides  of  a  stringer 
are  out  of  phase*  Zf  the  stringers  are  of  low  torsional  stiffness,  the  frequency  for 
the  panel  will  be  close  to  the  value  indicated  by  thin  Section  for  a  single  plate  with 
simply-supported  edges*  With  increasing  torsional  stiffness  the  lowest  resonant 
frequency  of  the  group  rises  to  approach  the  value  Indicated  by  this  Section  for  a 
single  plate  with  fixed  edges* 

The  highest  natural  frequency  of  a  panel  in  a  frequency  group  is  associated  with 
a  mode  of  vibration  in  which  the  half  waves  in  the  plates  on  the  two  sidos  of  a 
stringer  are  in  phase*  This  behaviour  corresponds  to  fixed-edge  conditions* 

4*4  Derivation 


4*4«Z  Timoshenko,  S. 
4*4*2  Warburton,  O.B. 
4*4*3  Clarkson,  B*L. 

4*4*4  Sewall,  J*L» 

4*^*5  Webster,  J.J. 

4*4*6  Petyt,  M. 

4*4*7  Webster,  J.J* 

Warburton,  G.B* 


Vibration  problems  in  engineering.  Second  edition. 

Von  Nostrand,  New  York,  1937* 

The  vibration  of  rectangular  plates* 

Proc.  I.Mech.E.,  Vol.168,  No. 12,  pp*37i-384,  1954* 

The  design  of  structures  to  resist  Jot  noise  fatigue* 

J.  R.  aeronaut.  Soc.,  Vol.66,  No. 622,  pp.603-6l6, 

October  1962* 

Vibration  analysis  of  cylindrically  curved  panels 
with  simply  supported  or  clamped  edges  and  comparison 
with  some  experiments. 

NASA  tech.  Note  D-379X,  January  1967* 

Free  vibrations  of  rectangular  curved  panels. 

Int*  J*  Mech.  Sci*,  Vol.10,  No. 7,  pp. 571-582,  July  1968* 

Vibration  of  curved  plates* 

J.  Sound  Vib.,  Vol.15,  No. 3,  PP*38l-395,  April  1971* 

Discussion  of  "Approximate  methods  for  the  determination 
of  the  natural  frequencies  of  stiffened  and  curved 
plates” • 

J.  Sound  Vib.,  Vol.18,  No.l,  pp. 139-141,  September  1971* 


4*5  Example 

It  is  required  to  estimate  the  lowest  natural  frequency  of  an  aluminium  alloy 
plate  in  both  the  simply-supported  and  fixed-edge  conditions.  The  plate  has  the 
following  dimensions  and  physical  properties » 


a  ■  300  mm, 

b  - 

200 

nun  9 

t  ■  0.8  mm, 

E  a 

70 

000  MN/m‘ 

Firstly 

a  222. 

b  "  200 

8 

1*5  * 

R  m  1000  mm, 
p  n  2770  kg/m^. 


2  2 
b_  200*  .n 

Rt  "  1000  x  0.6  *  30 


and 


V 


/ 


70  000  x  10u 
2770 


1 

50B0 


»  Q.99. 


4*5*1  All  edges  simply-supported 

Obtain  the  natural,  frequency  parameters  from  Figure  4*1* 


Mode 

a  n 

b2 

H* 

O 

• 

V*> 

2. 
n  t 

2*. 

f  »  VK  ~ 

m  b 

n*Rt 

m  »n 

b2 

m,n  ^2 

(m,n) 

(m/s) 

(nT1) 

(Hz) 

x.x 

1.3 

50 

12*9 

0.02 

255 

1,2 

3.0 

12.5 

2.87 

0.08 

227 

2, 1 

0*75 

50 

26*7 

0.02 

529 

1»3 

4*5 

5.536 

2.54 

0.18 

453 

3,1 

0*5 

50 

34.5 

0.02 

683 

28 


Hence  the  lowest  natural  frequency  with  all  edges  simply  supported  occurs  in  the  (l,2) 
mode  and  is  227  Hz. 

4 • 5 • 2  All  edges  fixed 

For  lowest  symmetric-symmetric  mode  from  Figure  4*2 

10"3  K  =  27.3  m/s  . 

For  lowest  symmetric-antisymmetric  mode  from  Figure  4.4 

10"3  K  »  20.0  m/e  . 

For  lowest  antisymmetric-symmetric  mode  from  Figure  4.6 

10“3  K  -  31.2  m/e  . 

For  lowest  antisymmetric-antisymmetric  mode  from  Figure  4*8 

10“3  K  ■  28.0  m/s  . 

Hence  the  lowest  natural  frequency  with  all  edges  fixed  occurs  in  the  first 
symmetric-antisymmetric  mode 

and  is  0.99  x  20.0  x  103  x  ■  396  Hz  . 

200 


KxlO* 

m/s 


FIGURE  4,8.  FREQUENCY  PARAMETERS  FOR  FIRST  ANTISYMMETRIC-ANTISYMMETRIC  MODE 
(PLATES  WITH  ALL  EDGES  FIXED) 


001 


FIGURE  4.12  SHAPES  FOR  SECOND  SYMMETRIC- ANTISYMMETRIC  MODE  (PLATES  WITH  ALL  EDGES  FIXED) 


FIGURE  4.14.  SHAPES  FOR  SECOND  ANTISYMMETRIC-SYMMETRIC  MODE  (PLATES  WITH  ALL  EDGES  FIXED) 


FIGURE  4  15.  SHAPES  FOR  FIRST  ANTISYMMETRIC-ANTISYMMETRIC  MODE  (PLATES  WITH  ALL  EDGES  FIXED) 
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Section  5 

THE  ESTIMATION  OF  R.M.S.  STRESS  IN 
STIFFENED  SKIN  PANELS  SUBJECTED  TO 
RANDOM  ACOUSTIC  LOADING 


) 

5.1  Notation  ‘ 


a 

width  of  skin  plate  or  frame  pitch 

a 

in 

b 

flat  panel  width  between  stringers  or 
curved  panel  arc  length  between  stringers 

a 

in 

E 

Young's  modulus  of  plate  material 

N/m2 

lbf/ in2 

f 

fundamental  natural  frequency  of  skin 
plate  assuming  all  edges  fixed 

Hz 

c/s 

G  (f) 
P' 

spectral  density  of  acoustic  pressures 
at  frequency  f 

(N/m2)2/Hz 

( lbf/ in2) 2/( c/i 

*5 

damping  ratio  correction  factor 

Lp.(f) 

spectrum  level  of  acoustic  pressure  at 
frequency  f 

dB 

dB 

P 

uniform  static  pressure  on  plate 

N/m2 

lbf/ in2 

^ros 

root  mean  square  fluctuating  pressure 

N/a2 

lbf/ in2 

R 

radius  of  curvature  of  panel 

a 

in 

So 

ratio  of  stress  at  rivet  line  in 
assumed  mode  shape  to  applied  uniform 
static  pressure  on  plate 

Srms 

r.m.s.  stress  at  rivet  line  due  to 
acoustic  loading 

N/m2 

lbf/ in2 

t 

thickness  of  panel 

n 

i*. 

V 

velocity  parameter  fur  plate  material4 

b 

damping  ratio  of  mode  of  vibration  of 
skin  panel 

P 

density  of  skin  material 

kg/m3 

ft 

<T  Poisson's  ratio  of  plate  material 

Both  SI  and  British  units  are  quoted  but  any  coherent  system  of  units  may  be 

used. 


3.2  General  Notes 


This  Section  gives  a  method  of  estimating  the  r.m.s.  stress  in  rectangular  skin 
panels  subject  to  random  acoustic  loading.  The  reference  stress  is  the  rivet  line 
stress;  stresses  at  other  positions  on  the  plats  may  be  found  by  applying  a  suitable 
factor  to  the  reference  stress. 

The  r.m.s.  streoo  for  a  stiffened  panel  subjected  to  random  acoustic  loading  on 
one  side  id  given  approximately  by  the  expression 

r  «  -.1/2 

Srae  "  [ft  f  °p<f)J  Vo*  ^ 


4  The  velocity  parameter  is  defined  in  SI  units  by  V  ■  (E/p)^/5080  or  in  British 

units  by  V  •  (E/p)^/200  000.  V  is  approximately  unity  for  all  common  structural 
metallic  materials. 

*  A  density  value  expressed  in  British  units  as  pounds  per  cubic  inoh  has  to  be 
divided  by  386.4  before  it  can  be  used  in  the  formula  for  V  given  here.  (A  force 

of  1  lbf  acting  on  a  mass  of  1  lb  produces  an  acceleration  of  386.4  in/s  •) 
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In  the  case  of  a  control  surface  or  similar  structure  where  two  skins  are 
coupled  by  ribs  and  both  vibrate  in  response  to  a  loading  largely  on  one  side*  the 
stress  level  is  about  l/3  of  the  value  given  by  Equation  (5*1) •  However*  where  the 
acoustic  loading  is  of  similar  intensity  on  both  aides  of  a  box  structure  the  stress 
level  is  about  l/2  of  the  value  given  by  Equation  (5«l)» 

In  deriving  Equation  (5«l)  it  has  been  assumed  that  the  predominant  form  of  skin 
vibration  is  one  in  which  individual  plates  within  the  stiffened  panel  vibrate 
Independently  in  the  assumed  mode  which*  for  this  calculation*  has  been  taken  to  be 
the  fundamental  fixed-edge  mode.  Thio  restricts  the  use  of  this  Section  to  skin- 
stringer  panels  where  the  bending  stiffness  of  the  stringers  is  sufficiently  high  to 
approximate  to  a  fixed-edge  condition  for  individual  plates.  It  is  not  applicable  to 
stiffened  panels  where  individual  plates  are  much  stiffer  and  ths  supporting  stiffeners 
relatively  more  flexible  than  for  conventional  skin-stringer  panels. 

The  pressure  is  assumed  to  be  uniform  and  in  phase  over  the  whole  of  each 
individual  plate  and  the  spectrum  level  of  the  acoustic  pressure  is  taken  to  be 
constant  over  the  range  of  frequencies  close  to  the  fundamental  natural  frequency  of 
the  panel. 

The  r.m.s.  stresses  given  are  the  nominal  unfactored  rivet  line  values  and  in 
using  these  stresses  to  calculate  a  fatigue  life  care  must  be  taken  to  ensure  due 
account  is  taken  of  any  geometric  stress  concentrations  in  the  region  of  the  rivet 
line. 


This  Section  may  be  used  for  etched  skins*  or  sK.ns  with  bonded  reinforcing*  by 
first  obtaining  the  nominal  rivet  line  r.m.s.  stress  for  the  panel*  assuming  it  to  be 
of  uniform  thickness*  and  then  applying  a  correction  with  the  aid  of  Reference  5*6'14. 

For  conventional  structures  without  special  damping  treatment*  f  may  be  taken 
as  the  undamped  natural  frequency. 


5 •  3  Notes  on  the  Figure-, 

Figures  5«1  and  5.2  givo  nomographs  for  S  .  The  nomograph  are  entered  at  a 

mo 

value  of  L  (f),  each  quadrant  being  used  in  turn  in  the  direction  indicated  through 
ps 

o 

ranges  of  a/b,  b/t,  b  /Rt  and  t/V.  Figure  5.2  is  a.-»  extension  of  the  stress 
range  of  Figure  5«1»  Figures  5.1  and  5*2  are  drawn  for  a  value  of  6  ■  0.017  ;  S  a 

for  other  values  of  6  is  calculated  using  correction  factor  which  is  plotted 

0 

against  6  in  Figure  5.3 


The  r.m.s.  stress  given  is  the  skin  surface  stress  across  the  rivet  line  at  the 
mid  position  of  the  longer  side  of  the  pla.e.  This  is  the  position  of  highest 
nominal  stress  when  the  plates  vibrate  in  their  fundamental  fixed-edge  mode.  No 
allowance  is  made  for  stress  concentration  factor  due  to  rivet  holes  etc.;  .his  is 
expected  to  be  negligible  at  distances  greater  than  2i5  times  the  rivet  di'^neter  from 
the  centre  line  of  the  rivet  along  the  rivet  line. 


In  the  derivation  of  the  nomographs  it  is  assumed  that  the  plate  bending  stress 
is  within  the  linear  region  where  it  is  directly  proportional  to  the  normal  prossure* 

that  is  p/E  less  than  about  20(t/b)**» 

The  higher  range  of  Figure  5*2  is  for  use  in  estimating  stresses  in  titanium  and 
high  strength  steel  panels.  However*  high  values  of  stresses  obtained  from  this 
Figure  should  be  treated  with  caution  as  no  extensive  measured  data  are  at  present 
available  to  check  the  validity  of  the  simple  response  theory  over  the  higher  stress 
range . 


5.4  Calculation  Procedure 


5.4.I  The  procedure  for  estimating  in  a  general  case  is  as  follows. 

nos 

(i)  Estimate  the  fundamental  natural  frequency  of  the  panel.  For  flat  plates 
use  Reference  5*6.13  or  Section  4  of  this  AGARDograph*  or  evaluate  the 
highest  frequency  of  the  first  group  of  natural  frequencies  using 
Section  3  of  this  AGARDograph.  For  curved  plates  calculate  the  natural 
frequency  using  Section  4  of  this  AGARDograph* 

(ii)  Obtain  the  value  of  spectrum  level  of  acoustic  pressure  L  (f)  at  the 

ps 

calculated  frequency.  If  only  the  band  pressure  level  is  known*  it  is 
first  corrected  to  pressure  spectrum  level  (unit  bandwidth)  using 
Reference  5*6*11*  The  reference  pressure  for  sound  pressure  level  is 

20pN/m2. 
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(iii)  Calculate  the  parameters  a/b,  b/t>  b  /Rt  and  t/V  and  read  the 
value  of  S  from  Figures  5*1  or  5*2. 

XTOC 

(iv)  For  values  of  6  other  than  0.017»  factor  the  estimated  value  of 

S  by  K.  obtained  from  Figure  5«3*  The  value  of  6  =  0.017  i® 
rrns  o 

typical  of  aircraft  structures  without  special  damping  treatment. 

5« 4*2  Within  the  nomograph  the  spectrum  sound  pressure  level  is  converted  into 
the  spectral  density  of  acoustic  pressure.  The  spectrum  sound  pressure  level  is 

converted  into  the  root  mean  square  fluctuating  pressure  in  units  of  (N/m  )/Hz  (see 
expression  below  or  Reference  5.6.12)  and  then  squared  giving  a  value  in  units  of 

(N/m  )  /Hz*’.  Since  unit  bandwidth  is  used  this  is  numerically  equal  to  the  spectral 

density  of  acoustic  pressure  Gp(f)  in  units  of  (N/m  )  /Hz. 

Lp8(f)  -  20(log10  Prm8  ♦  4.70)  . 

If  L  (f)  is  required  in  British  units  of  (lbf/in^)^/(c/a)  it  is  given  by 

P  • 

Lps(f)  "  20( log^Q  Prm8  ♦  8.54)  • 

5*3  Comparison  with  Measured  Data 

Figure  5«4  shows  a  comparison  of  estimated  and  measured  stresses  in  flat  or 
curved  plates.  Figure  5*5.  shows  a  comparison  of  estimated  and  measured  stresses  in 
control  surfaces  and  other  coupled  skin  box  type  structures.  For  these  estimated 
values  the  factors  given  in  Section  5*2  were  applied  to  the  stresses  given  by  Equation 
(5«l).  In  estimating  r.m.s.  stresses  to  compare  with  measured  values,  6  has  been 
assumed  to  be  0.017. 
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5.6.13  -  Natural  frequencies  of  uniform  flat  plates. 

Engineering  Sciences  Data  Item  No.  66019, 
February  1966. 

5.6.14  -  The  effect  of  edgo  reinforcement  on  the  stresses 

in  skin  panels  undor  uniform  pressure. 
Engineering  Sciences  Data  Item  No.  67029. 

May  1967. 

3*7  Examples 


5> 7*1  It  is  required  to  estimate  the  rivet  line  r.m.s.  stress  in  a  stiffened 
panel  subjected  to  Jet  noise  on  one  side.  The  variation  of  found  pressure  level  over 
a  range  of  frequencies  is  given  in  the  table,  sound  pressure  levels  being  l/3  octave 
band  levels. 


Sound  pressure 
lovel (dB) 

130 

145 

150 

148 

143 

130 

Frequency  (Hz) 

150 

200 

300 

500 

700 

1000 

The  panel  is  made  up  from  uniform  plates  having  the  following  dimensions  and 
properties  1 

a  »  210  mm,  b  »  140  mm,  t  «  1.2  mm,  R  ■  1500  mm, 

E  »  70  000  MN/m2,  p  »  2770  kg/m3 * 5,  5  «  0.020  . 


Firstly 


a 

b 

b 

t 


210 

140 


1.5 


-  116.7  , 


£_ 

Rt 


140 


1500  x  1.2 


10.9 


70  000  x  10 
2770 


x  5^50 


0.99 


and 


t 

V 


1.2 

0.99 


a  1.21  mm. 


From  Section  4  of  this  AGARDograph  the  fundamental  natural  frequency  of  the 
plates  with  all  edges  fixed  is  obtained. 

f  «  0.99  x  10  200  x  -  618  Hz. 

140 

By  interpolation  from  the  table  the  1/3  octave  band  pressure  lovel  at  618  Hz  is 
145*7  dB. 

From  Reference  5*6.11 

L  (f)  »  145*7  -  21.7  -  124  dB. 

ps 

2 

From  Figure  5*1.  entering  the  nomograph  at  124  dB,  S__.  »  14*0  MN/m  for 

rnj" 

5  «  0*0x7* 

From  Figure  5*3  for  6  *  0*020f  «  0*922*  Hence  the  rivet  line  r*m*0* 

stress  =  14*0  x  0.922  «  12.9  MN/m2. 


5*7*2  If  a  doubler  is  now  bonded  to  the  skin  at  the  rivet  line  the  new  stresses 
may  be  found  usinG  Reference  5*6*14*  Additional  definitions  in  this  Reference  are* 


br 

width  of  reinforcement 

m 

in 

fo 

nominal  stress  at  rivet  lino  in 
unreinforced  panel 

N/m2 

lbf /  in' 

fl 

nominal  stress  at  rivet  line  in 
reinforced  panel 

N/m2 

lbf/in' 

f2 

nominal  stress  at  edge  of 
reinforcement 

N/m2 

lbf/ in‘ 

*r 

total  thickness  of  panel  and 
reinforcement 

m 

in 

In  this  example 

f  is  the  previously  calculated  value  of 

S 

rmo 

The  doubler  in  the  same  material  as  the  plate  and  has  the  following  dimensions t 

b  ■  21  mm  and  t„  »  1.8  mm* 


r 

r 

Firstly 

b 

r 

b 

21 

*  140  “ 

0.15 

and 

*r 

t 

1.8 

*  1.2  “ 

1*5* 

From  Reference  5*6.14 

fl 

fo 

a  0.485 

and 

—  «  0.677. 

o 

Hence  the  rivet  line  r.m.s.  stress  with  the  doubler 

«  12*9  x  0.485  -  6.3  MW/m2 

and  the  r.m.s.  stress  at  the  edge  of  the  doubler 

»  12.9  x  0.677  =>  8.7  MN/m2. 

It  should  be  noted  that  Reference  5*6.14  is  strictly  applicable  to  panels  of 
a/b^2.  Within  the  accuracy  of  the  simple  theory  for  stress  response  to  acoustic 
loading  it  may  bo  used  for  lower  values  of  n/b  . 
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fef/in* 


10 


I  4  91  I 


3  4  3  «  t 


10 


Syn*  MN/m* 


yo  dB 

FIGURE  5.2.  HIGH  RANGE  STRESS  NOMOGRAPH 


FIGURE  5.3  OAMPING  RATIO  CORRECTION  FACTOR 
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Estxnoted 

MN/m* 


I  10  itf 


Mflosured  Sf!nj  MN/m* 


FIGURE  5.4.  COMPARISON  OF  ESTIMATED  AND  MEASURED  STRESSES 
IN  FUAT  AND  CURVED  PANELS 


